Composite Geometric Phase for Multipartite Entangled States 
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When an entangled state evolves under local unitaries, the entanglement in the state remains 
fixed. Here we show the dynamical phase acquired by an entangled state in such a scenario can 
always be understood as the sum of the dynamical phases of its subsystems. In contrast, the 
equivalent statement for the geometric phase is not generally true unless the state is separable. For 
an entangled state an additional term is present, the mutual geometric phase, that measures the 
change the additional correlations present in the entangled state make to the geometry of the state 
space. For A*' qubit states we find this change can be explained solely by classical correlations for 
states with a Schmidt decomposition and solely by quantum correlations for W states. 
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In a seminal paper [l|, Berry recognized that a quan- 
tum system undergoing a cyclic, adiabatic evolution 
records the path of its evolution as a geometrical quan- 
tity in the phase of its wavefunction. This phase, the 
geometric phase (GP), forms part of the total phase of 
wavefunction in addition to the more familiar dynami- 
cal phase. Since Berry's initial discovery, the GP has 
been found to occur in more general circumstances; non- 
adiabatic and non-cyclic evolutions 0], non-abelian 
form Q and for states that are mixed 

GPs are of interest for many reasons, amongst them 
topological effects in many body systems ^ and their 
use for quantum information processing JTj, a paradigm 
where entanglement is known to be key in its advantage 
over the classical counterpart. GPs may also be used 
as a tool. Because they depend on the geometry of the 
space the states traverse, GPs can provide information 
on this space. As an example, this property can be used 
to discover the coordinates of quantum phase transitions 
in a parameter space 0]. On another currently popu- 
lar front, entanglement has been shown to be present in 
many body systems at finite temperatures ^9]. Studying 
GPs of multipartite entangled states may therefore prove 
to be complimentary to both efforts. 

GPs have been studied only for entangled bipartite sys- 
tems so far. These include qubits precessing in magnetic 
fields Td^l and general evolution jT3| both without inter- 
action (fixed entanglement) and various specific Hamil- 
tonians for bipartite systems with interactions (changing 
entanglement) [H, [IjI- It has been shown [l^, |Tl[ in 
bipartite systems that even if there are no interactions 
during the evolution, fixed entanglement affects the geo- 
metric phase. 

It is this last relation that we investigate further in this 
Letter for multipartite systems. Of the two terms mak- 
ing up the total phase of a quantum state, dynamical and 
geometric phase, we will show that under local evolutions 
the dynamical phase of a composite state can always be 
understood as the sum of the dynamical phases that com- 
prise the composite system whether the state is entangled 



or not. We will also show that the same statement can- 
not be made for GP when the state is entangled. When 
the composite state is entangled an additional correction 
term, much like an interference term, A7, is present. Ex- 
plicitly r = Y^n=i + ^7 where the composite state 
GP is r and the n*'' of the N subsystems has GP 7,^^. 
The correction term we call the mutual GP. Each sub- 
system, pm is obtained by tracing over the other — 1 
subsystems in the composite state p. 

In tracing out the subsystems from the composite state 
we have removed both quantum and classical correlations 
14| . The extra correlations present in the composite 
state, unaccessible when one only has access or control 
on one of the subsystems, modify the GP but not the 
dynamical phase when the state is entangled. Should 
one only have access or control on one of the subsystems, 
Pn, then a GP of 7^^ would be observed. As previously 
noted GPs are dependent on the underlying geometry of 
the quantum state space and here we are using the GP as 
a tool to sample the change these correlations have on the 
geometry of the state space. We characterize which type 
of correlations, quantum (entanglement) or classical, are 
responsible for the correction, A7, for two specific types 
of state we know how to quantify the entanglement of 
using the relative entropy of entanglement. En [l^. We 
will show A7 can be attributed solely to classical corre- 
lations for states with a Schmidt decomposition (labeled 
hereafter as S states) and solely to entanglement for W 
states. 

Our aim is to study the effect of fixed entanglement 
on the mutual GP, A7. We fix entanglement during the 
evolution by requiring that the evolution on the compos- 
ite state, the unitary lA, be composed of local unitaries, 
Un, acting on each subsystem, n. We write this as 



N 



U{t) 



I Unit). 



(1) 



n=l 



With this condition we can show that the dynamical 
phase of the composite system. A, is always given by 
the sum of the dynamical phases of its subsystems, . 
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The definition of dynamical phase [1] is 

A = -i tr[pU\t)U{t)]dt. 



(2) 



Substituting eq ([T]) into this definition we see 



(3) 



showing A — J2n=i ■ The composite dynamical phase 
is always equal to the sum of the subsystem dynamical 
phases whether the state is entangled or not. This state- 
ment holds because we have constrained the dynamics 
to be local. Before demonstrating this is not the case 
for GP we discuss parallel transport as a useful way to 
isolate the GP. 

A quantum state is said to be parallel transported 
when it acquires no dynamical phase at each point along 
its evolution. Formally the mathematical condition that 
requires the state is in phase with itself at each point is 
trpU{t)^U{t) = Oy t. If the state is parallel transported 
(by W'l ) then the total phase obtained by the state will be 
equal to the GP. Pancharatnam gave a natural pre- 
scription to obtain the total phase, generalized to mixed 
states [1] . As mentioned, under parallel transport condi- 
tions this becomes the GP, F, accumulated over the time 

te[0,T]. 



F = argitrpUHr)}. 



(4) 



For unentangled states of the form p = ^^=1 

see that under local unitaries eq ([4]) becomes F = 

arg{nlii'^PnC/|(r)} so F = E1i7^'- But in gen- 

eral we see that F / J2n=i ■ 

It is known that there are many parallel transport con- 
ditions for mixed states. In this study the subsystems, 
p„, will be the mixed states. Some of these parallel trans- 
port conditions produce GPs that are a property of an 
arbitrary entangled purification of the mixed state in a 
higher dimensional space, not just of the mixed state it- 
self. It is also known, however, that a subset of these con- 
ditions, those proposed by Q ('stronger' parallel trans- 
port conditions) provide a mixed state GP that is a prop- 
erty of the evolution of the mixed state only 17| . For this 
reason we use the stronger parallel transport conditions 
to constrain the local unitaries. These conditions require 
that every eigenvector, of each subsystem, n, be 

parallel transported, formally |[/t(i)?7„(t)|0^) = V 
i,n,t. We write the unitary acting on subsystem n that 



fulfills these conditions U, 



We now calculate F and 7^^ for arbitrary superposi- 
tions of W states which contain, amongst others, W, S 
and GHZ states. We will write these states as 



N 



|vI'(0))=^afe|iV,fc), 



(5) 



where ak are the complex probability amplitudes and the 
W state \N, k) is defined by 



|iV,fc) 



1 
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(6) 



S is the total symmetrization operator. Time evolution 
of the state |*(0)} is given by |*(i)} = U{t)\^{Q)). Each 
individual qubit in the state considered on its own is 
given by the density matrix /9„ (i) = tri,,N^^n\^{t)){^{t)\- 
In the |0), |1) basis the subsystem state explicitly is 

P„(0)=poo|0)(0|-f pn|l)(l|+poi|0)(l|+pSi|l)(0|, (7) 



Pii - Etok^Pw and 



N ' 



where poo - 

Poi = ELi al.^au^KN ~k + 1)/N. The GP for 
both the composite state \'^{T)) and each of its sub- 
systems, pn{T), is made by substitution into eq 
We have F = a^Tg{J2k,i=oa.iak{N,l\uHT)\N,k)}. Cross 
terms disappear from this equation when each local uni- 
tary, Un, brings each subsystem back to the same ray 
(C/|(T)|0} = C/|(0)|0) up to a global phase). Under this 
evolution the global GP, F, becomes 



arg • 



N 



k=0 ( 



(8) 



A:=0 



expressing the composite system GP in terms of pure 
state qubit phases, 7„. We have used the polar represen- 
tation to define c„e''''" := (0|t/„(T)|0) and its complex 
conjugate c„e~*^" := (l|[/„(r)|l). 7„ (-7„) arc the GP 
that the pure qubit states |0) (|1)) acquire over the evo- 
lution Un(T) and Cn = 1 when each subsystem unitary 
is cyclic. The matrix A'' has elements 1(— 1) to capture 
the sign of 7„ (— 7ri) in the exponent. Each ordered row 
in A'' has N — k elements with value 1 and k elements 
of —1. There are ( " ) ordered rows in A^ given by all 
permutations of the first row elements. 

Eq ([5]) can be understood as follows: Each W state, 
|iV, fc), is a superposition of ( " ) kets, each with phase 

factor e'^"=i The phase factors for each ket are 

averaged in the m summation resulting in an overall 
phase factor for each W state (fc). The phase factors 
for each W state are then averaged (the k summation) to 
give the total overall GP factor of the state. 

The states of the subsystems are given by p„ and the 
GP associated to the them is also obtained by substi- 
tution into eq ([4]) with [/|(T). We substitute the gen- 
eral form of the eigenvectors and eigenvalues in to this 
equation to obtain 7,f = arg{i(l -I- r){(j)1\ul{T)\(j)'l) + 
i(l — r){(j)2\u\{T)\(j)2)} , where the general eigenvectors 
of a two level system are |(^") = e^^^ cos(0„/2)|O) -|- 
sin(6'„/2)|l) and |<^5) = -e^ sin(0„/2)|O) + 
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cos(0n/2)|l). Rewriting in the |0), |1) basis and 
substitution of 7„ identities results in the cychc GPs of 
the subsystems to be 7*^ = arg{cos7„ + ircos0„ sin7„}. 
This can also be written as 



7*^ — arctan 




, N -2k 
N 



tan 7„ 



(9) 



for non-degenerate eigenvalues of /9„. We can in- 
terpret this result in a similar manner to the 
pure state result. Writing eq (O as 7*^ = 



arg 



2 Af-fc z-Y„ 
1 N ^ 



see 7^ is an average of the phase factors of the states 
0) and |1) weighted by their relative frequency in the 
composite state. For the remainder of this Letter we will 
discuss specific cases of eqs ([8]) and dH). 

To simplify things, each local unitary will now be iden- 
tical i.e. Wll = C/ll in aU the following. All pure 
state qubit GPs, 7„, also become identical and will be 
labeled 7. 

The first specific example is that of S states. For these 
states only two amplitudes in eq ([5]) are non-zero, oq and 
UN. The S state is 



l*(0)) - yi(l+r)|0)®^ + yi(l-r)|l)«^. (10) 

We have written oq and cat as the square root of the 
eigenvalues of a bipartite state. Substitution into eq ([8|) 
yields 

r= arg{i(l-|-r)e'^''-t-i(l-r)e-*^^} 
= arctanjrtan A^7}. (11) 

In this formula pure state phases rather than phase fac- 
tors add up. Each of the local GPs for the S state are 



7 



M 



arctan{rtan7}. 



(12) 



To quantify the amount of entanglement in the composite 
states we use the relative entropy of entanglement, Eji 
Hill. For S states Er is given by 



En = l--[il + r) log2 (1 + r) + (1 - r) log^ (1 - r)] . (13) 

Efi takes a maximum value of 1 when r = and a mini- 
mum value when r = 1. When a S state is maximally 
entangled (r = 0) F can only take the values and tt, 
however eq ([12]) for j'^ is no longer valid as p„ is degen- 
erate and must be evaluated via path ordering (see [l9| ) . 
At the other extreme, separable states (r ~ 1), we see 
F = as expected. 

The second specific example is a W state |^'(0)) = 
|iV, k). The identical local evolution in this case is par- 
ticularly appropriate as these states often occur when 



subsystems are indistinguishable. For this state |afcp = 1 
and the composite GP takes a particularly simple form 



F = (iV - 2fc)7. 



(14) 



For a W state the GP is the sum of the pure state qubit 
phases much like a pure separable state. This result in 
fact holds for any arbitrary superposition of W kets (for 
example ai|001)+a2|010)-|-a3|100}). However the subsys- 
tems are not generally pure and their GPs are generally 
not 7 and —7. They are 



A f 

7 = arctan 



A^-2fc 
TV 



tan 7 



(15) 



The relative entropy of entanglement for W states is 

Mm 



known to be 



Er. 



log2 



^(7V-fc)l0g2 



N-k \ (k 
(16) 

Entanglement is maximal when iV = 2fc in which case 
Eji — N — log2 {1^/2) and is minimal {E^i — 0) when 
k — 0,N (the state is separable). When the W state is 
maximally entangled F = 0. This limit is similar to the 
maximally entangled S state, except it could also take the 
value TT. We also note, for the same reason as maximally 
entangled S states, eq (fTS]) is not valid for calculation of 
7^^ when N — 2k. Again F = A^7*^ for separable states 
(when k^O,N). 

We now have the ingredients to calculate mutual GP, 
A7, written explicitly as 



N 



A7 



n=l 



M 



(17) 



As previously mentioned, we have removed the quantum 
(entanglement) but also the classical correlations in trac- 
ing out each subsystem from the composite. We can char- 
acterize which missing correlations are responsible for the 
correction term A7 using ideas from entanglement dis- 
tance measures, specifically here, the relative entropy of 
entanglement, Efj. 

To calculate Efi, the closest separable state, ps, to the 
entangled state needs to be found. Here closest means 
the separable state from the set of all separable states 
that minimizes the quantum relative entropy S{p\\ps) — 
tr{plog p — plogps), between it and the entangled state, 
p. The minimum value of S{p\\ps) is equal to Efi. The 
procedure maximizes classical correlations between the 
entangled and the separable state, any correlations left 
over are quantum. The closest separable states are given 
in [lil for S states and ^ ^ for W states. The GPs, 
F5, for these closest separable states are 

Ts = arg{l(l + r)e»^^ + 1(1 - r)e-'^n (18) 
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for the S state and 

Ts = TVarctan | tan7| (19) 

for the W state. In the latter we have used the binomial 
theorem and the property N aig{a} = arg{a^}. 

By taking the difference between the GPs of the com- 
posite state and the closest separable state (F — Fs) we 
can see exclusively what difference the quantum correla- 
tions make to the GP. For S states this difference is zero 
as F = Fs- We can therefore state quantum correlations 
have no effect on GP and state space geometry. We can 
also state the correction term A7 is the sole result of clas- 
sical correlations for S states. In contrast, F — Fg = A7 
for W states because Fg = J2n=i ■ Here we can state 
A7 results solely from quantum correlations. One com- 
mon feature the states share is that entanglement does 
not affect curvature when Er < 1. Note however that 
S states can have a continuum of values < E^j < 1 
independent of N. E^ of W states does have N depen- 
dance but discrete values. The smallest of these values 
are Efj = (separable, k = 0,N), E^ = 1 (singlet state, 
N — 2,k = 1) or Er > log2 e for fc = 1, iV— 1. According 
to the relative entropy of entanglement, perhaps counter 
intuitively, W states are more quantum than S states, in 
the sense that classical correlations cannot approximate 
these states as closely. This is related to the robustness 
of the entanglement. Almost all correlations in S states 
are classical, tracing out just one subsystem results in a 
separable state. In contrast W states remain entangled 
to the last pair of subsystems. Plotted in figure [T] is A7 
for S and W states, the difference classical and quantum 
correlations make to the GP respectively. Although we 
don't know how to quantify the entanglement of states 
like eq (O we expect that A7 for these states can be 
attributed to a mixture of classical and quantum correla- 
tions, each reducing to what are probably the two special 
limiting cases, S and W states. 

In this Letter we have demonstrated that under lo- 
cal evolutions, entangled states gain a correction term 
to their GP not present in the dynamical phase. This 
correction term A7 has its origins in the change the ex- 
tra correlations present in entangled states make to the 
state space geometry and we have showed it can solely be 
attributed to classical correlations for S states and quan- 
tum correlations for W states. Here we have used the 
relative entropy of entanglement definition of quantum 
and classical correlations and it would be interesting to 
see what statements can be made for other entanglement 
measures. GPs are path dependent quantities but by de- 
termining this path dependence one can learn something 
of the geometry of the underlying space. We hope that 
GPs may prove useful in determining how classical and 
quantum correlations modify this geometry. 

MSW & VV acknowledge funding from EPSRC, 
QIPIRC and the Royal Society. 
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FIG. 1: A7 plotted against entanglement, Er, and 7 for 
A*' = 51. Left panel: S state (classical correlations only). 
Right panel: W state (quantum correlations only). A7 is 
non-modular and plotted in the large limit with varying 
path, 7, and entanglement in the composite state, Er. For 
constant 7 S states are monotonic w.r.t. Er. W states are 
monotonia except for the region \'y\/iT < 1/2. Periodicities 
w.r.t. 7 are n and 2ir for S and W states respectively when 
Er is fixed. For both states A'^ controls the magnitude of A7. 
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